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Systematic studies’ * have revealed hundreds of ultra-compact dwarf galaxies (UCDs°) 
inthe nearby Universe. With half-light radii r, of approximately 10-100 parsecs and 
stellar masses M,~ 10°-10° solar masses, UCDs are among the densest known stellar 
systems. Although similar in appearance to massive globular clusters’, the detection 
of extended stellar envelopes**’, complex star formation histories”, elevated 
mass-to-light ratio™” and supermassive black holes ™ suggest that some UCDs are 
remnant nuclear star clusters” of tidally stripped dwarf galaxies”, or even ancient 
compact galaxies”. However, only a few objects have been found in the transient 
stage of tidal stripping”, and this assumed evolutionary path” has never been fully 
traced by observations. Here we show that 106 galaxies in the Virgo cluster have 
morphologies that are intermediate between normal, nucleated dwarf galaxies and 
single-component UCDs, revealing a continuum that fully maps this morphological 
transition and fills the ‘size gap’ between star clusters and galaxies. Their spatial 
distribution and redder colour are also consistent with stripped satellite galaxies on 
their first few pericentric passages around massive galaxies”. The ‘ultra-diffuse’ tidal 
features around several of these galaxies directly show how UCDs are forming through 
tidal stripping and that this evolutionary path can include an early phase as a nucleated 
ultra-diffuse galaxy”. These UCDs represent substantial visible fossil remnants of 
ancient dwarf galaxies in galaxy clusters, and more low-mass remnants probably 
remain to be found. 


We performed a detailed structural analysis (Methods) of compact 
stellar systems in the nearby Virgo cluster of galaxies based onthe Next 
Generation Virgo Cluster Survey (NGVS). Two intermediate populations 
were identified between typical nucleated dwarf galaxies (dE,Ns) and 
ultra-compact dwarf galaxies (UCDs). First, we have selected 51dE,Ns 
that have ‘overweight’ central nuclear star clusters (‘strongly nucleated’ 
dEs), with a nuclear-to-total-luminosity fraction of fauc > 8%, that are 
among around the highest 10% in fauc for about 600 Virgo dE,Ns with 
g’-band absolute magnitude M, between -16 and -9.5 mag. Second, we 
found that approximately 15% of Virgo UCDs* cannot be described by 
a single King or Sérsic profile with reasonable parameters, indicating 


the existence of extended envelopes. Follow-up Gemini Multi-Object 
Spectrograph (GMOS) spectroscopy (Methods), together with previ- 
ous observations”©””, leads to a clean sample of 55 UCDs with envelopes 
(eUCDs) in Virgo, of which 50 were selected with NGVS data (including a 
few bright eUCDs identified in previous Hubble Space Telescope (HST) 
work®"*"5) and an additional five eUCDs that could only be resolved in 
the HST imaging from the Advanced Camera for Surveys Virgo Cluster 
Survey (ACSVCS) (Methods). Our sample fully bridges the morphologi- 
cal gap between dwarf galaxies and UCDs, as shown in Fig. 1a, despite 
the fact that eUCDs are much more compact than strongly nucleated 
dEs and were not initially discovered as galaxies. 
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Fig. 1| Strongly nucleated dEs, eUCDs and UCDs in the Virgo Cluster core. 
a,u`g'i’' colour composite images of representative strongly nucleated dEs (left 
column) and eUCDs (right column). Each cutout image is 0.5 arcmin x 0.5 arcmin. 
b, BSDVS imaging of low-surface-brightness features in the Virgo core region. 
Yellow and dark red circles represent UCDs and eUCDs, respectively. Blue circles 


Are these strongly nucleated dEs and eUCDs truly transition 
objects? The imaging depth (surface brightness 4y ~ 29 mag arcsec’) 
of both NGVS and the Burrell Schmidt Deep Virgo Survey (BSDVS)?8 
has enabled us to visually identify notable S-shaped tidal features 
and streams—unique indicators of the ongoing tidal disruption” — 
around a few strongly nucleated dEs and eUCDs (Fig. 1b, see Meth- 
ods for more details). Only a few Virgo galaxies with similar features 
have been reported before”!”5“°, We find that VCC 1672 near M89 and 
NGVS J123037.24 + 124609.2 (hereafter NGVS 2078) near M87 both 
show tidal features at 4y > 27-28 mag arcsec” that span more than 
3 arcmin (about 14 kpc). We also find the approximately 20-kpc-long 
tidal streams associated with eUCD NGVS-UCD156 (near NGC 4365) 
and the bright globular cluster NGVSJ122846.19 + 134311.4 (hereafter 
NGVS J1229 + 1343)*, and a tentative 70 kpc stream originating from 
eUCD NGVS-UCDS09 near M87. Locations and zoom-in images of tidally 
disrupting galaxies in the Virgo core region are highlighted in Fig. 1b. 
Their proximity to massive galaxies and resemblance to the distorted 
morphologies shown in numerical simulations”*”’ strongly suggest 
ongoing tidal disruption rather than imaging artefacts. These new 
findings bring the number of known disrupting dE,Ns in Virgo into the 
range predicted by simulations”. 

On the size-luminosity scaling relationship (Fig. 2), most strongly 
nucleated dEs with M, < -11 are actually larger and thus more diffuse in 
surface brightnesses relative to the main locus of Virgo dE,Ns, suggest- 
ing tidal interaction and heating to puff up galaxies. Only the fainter 
(Mg = -10) strongly nucleated dEs have sizes consistent with their 
normal dE,N counterparts. Figure 2 also shows that the sizes of eUCD 
envelopes span approximately 30-200 pc, filling the long-observed 
size gap between star clusters and galaxies”. Strongly nucleated dEs 
and eUCD envelopes map out a contiguous sequence in decreas- 
ing size and luminosity, something which is expected by repeated 


and squares are strongly nucleated dEs with and without velocity measurements. 
u’g’i’ co-added and Gaussian-smoothed deep NGVS imaging of four disrupting 
dE,Ns (VCC 1672, NGVSJ 1229 + 1343, NGVS 2078, VLSB-A) and NGVS-UCDS09 
are shown in the insets. The arrow indicates the 70-kpc-long tidal stream that 
probably originates from NGVS-UCDS509. 


episodes of tidal stripping at each pericentric passage of a massive 
host”. 

Those dE,Ns undergoing active tidal disruption (VCC 1672, NGVS 
2078, NGVS J1229 + 1343 and VLSB-A) show more complicated struc- 
tures, which includes a compact inner stellar envelope with effective 
radius R, of a few hundred parsecs surrounding the central nuclear 
star cluster, and an elongated and twisted outer stellar component 
or stream with (He » > 27 mag arcsec’. Such ‘tidal breaks’ in surface 
brightness profiles, which has been predicted in simulations”, requires 
athree-component modelto describe (Methods), rather than the clas- 
sical two-component ‘nucleus + main galaxy’ model for undisturbed 
dE,Ns. The size-luminosity diagram shows that the main galaxy 
(inner + outer) components (pentagons in Fig. 2) occupy the same 
parameter space as Virgo ultra-diffuse galaxies (UDGs)*°, while the 
central nuclear star clusters (NSCs) (circles in Fig. 2) are indistinguish- 
able from UCDs and massive globular clusters (GCs). However, if the 
outer component is entirely removed (or is below the surface bright- 
ness detection limit), the inner galaxy component alone (diamonds 
in Fig. 2) is structurally similar to strongly nucleated dEs or eUCDs. 
These components are connected in Fig. 2 and provide details about the 
potential evolution in the size-luminosity diagram, ending up as UCDs. 

Previous simulations®?** show that the change in stellar mass and 
half-light radius of a tidally stripped galaxy can be parameterized as a 
function of the total amount of mass lost within the half-light radius and 
the dark matter halo density profile. Figure 2 plots these evolutionary 
tracks of total stellar mass and half-light radius for a typical satellite 
(R= 1kpc, M, ~ 108° M,) that is embedded in a cored (upper track) or 
cuspy (lower track) dark matter halo profile. Ifthe progenitors were dEs, 
as proposed here, cored dark matter halos may more easily produce an 
evolutionary stage with large and diffuse envelopes, consistent with 
UDGs*. More extended disky progenitors would be needed to produce 
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Fig. 2 | The size-luminosity scaling relationship of Virgo early-type stellar 
systems. The stellar envelopes of strongly nucleated dEs (blue) and eUCDs 
(dark red) form a distinct population when compared to other Virgo nucleated 
galaxies (grey diamonds) from NGVS*, with eUCDs filling the size gap between 
normal galaxies and star clusters**”*°. Some strongly nucleated dEs have sizes 
comparable to UDGs (dark/light green for non-nucleated/nucleated UDGs*”). 
Tidally disrupting dE,Ns are decomposed and each labelled with three 


sufficient UDGs in the case of cuspy dark matter profiles”. The diversity 
in size and surface brightness of strongly nucleated dEs suggests that 
their dark matter halo profiles might differ. Those dE,Ns possessing very 
diffuse stellar envelopes that are comparable to ultra-diffuse galaxies 
(UDGs”) are more likely to be embedded in cored dark matter halos”, 
while cuspy dark matter halos can account for those more compact 
and higher surface brightness envelopes. The existence of very diffuse 
disrupting dE,Ns provides more direct evidence that tidal stripping 
contributes to the production of both the UCD and UDG populations”. 
The evolutionary pathway to forming UCDs may sometimes include 
a transient stage as a nucleated UDG, which through stripping would 
also be dark matter deficient”. 

Tidal stripping drives the increase of the nuclear-to-total-luminosity 
fraction (fruc), as NSCs are much more resilient to tidal stripping and will 
remain mostly unchanged until the entire galaxy is fully disrupted”. 
For atypical dE,N ( fauc = 1%) to evolve into one of our strongly nucleated 
dEs (foue ~ 8-35%) via tidal stripping, it would lose more than 90% of 
its stellar mass. To form a eUCD with even higher f uc (about 30-80%), 
it will have lost about 99% of its stellar mass. For the average eUCD 
((2’ toeucn = 20.0 mag and ¢ fnucdeucn = 56%), this corresponds to atypical 
dE,N progenitor of approximately 10° M,. 

Compact stellar systems are predominantly more metal-rich and 
redder than typical galaxies at the same stellar mass”. Figure 3 shows 
the colour-magnitude diagram for stellar envelopes of eEUCDs and 
strongly nucleated dEs, compared to the inner colour (within effective 
radius) of Virgo early-type galaxies*®. The mean envelope (g’-i’) colours 
of eUCDs and strongly nucleated dEs are respectively about 0.07 and 
0.03 mag redder than galaxies at the same magnitude, which is similar 
to the colours of dE,Ns with 17 < g’ < 19 mag (M, = 10”° M,), 2-4 mag 
brighter than the bulk of eUCDs. This provides independent support 
for the idea that these eUCDs started out as galaxies approximately 
10-40x more massive, with more normal f,,,.. The smaller red offset for 
the strongly nucleated dEs is consistent with their lower fauc both of 
which indicate that these galaxies have not undergone as much trans- 
formation as the eUCDs. A fully stripped dE,N will contribute its NSC 
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connected points representing each component. The evolutionary tidal 
stripping track™ for galaxies in cored (upper) and cuspy (lower) dark matter 
halos are shownas dark blue curves with arrows denoting where the track is 
extrapolated. The remaining mass fraction, compared to the initial stellar mass 
(m/m,), is labelled along the track. Dashed lines from lower-left to upper-right 
denote constant mean effective surface brightness. A representative loerror 
bar is shown inthe legend. 


to the large population of UCDs (g’ = 20-21 mag) or globular clusters 
(g’ = 22-23 mag or fainter). We found that the nuclear colour of eUCDs 
and bright strongly nucleated dEs shows no difference relative to the 
colour-magnitude relation of galaxy NSCs“ (Extended Data Fig. 7). 
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Fig. 3 | Colour-magnitude diagram for eUCDs, strongly nucleated dEs and 
Virgo early-type galaxies. The stellar envelopes of eUCDs and strongly 
nucleated dEs are redder than their dE counterparts at the same luminosity. 
The mean envelope (g’-i’) colour of eUCDs (red horizontal band) is about 

0.07 mag redder and the mean envelope (g’-i’) colour of strongly nucleated dEs 
(blue horizontal band) is about 0.03 mag redder than the dwarf galaxy colour- 
magnitude relationship (dark blue line) at the same magnitudes. The thicknesses 
of the red and blue bands represent the louncertainties in the mean envelope 
colour for eUCDs and strongly nucleated dEs, respectively. Where these bands 
intersect the colour-magnitude relationship shows what the mean progenitor 
luminosity for each population is expected to be, assuming evolution due to tidal 
stripping. A loerror bar for the colour and magnitude is shown for each object. 
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Fig. 4 | Morphological sequence of tidal stripping. Each 3.75 arcmin x 3.75 

arcmin panel shows adE,N or a UCD, ordered by increasing nuclear luminosity 
fraction. The sequence starts with VCC 1107, which is an undisturbed dE,N at the 
starting point of the evolutionary track in Fig. 2, and the subsequent galaxies all 


The structural properties and colours of strongly nucleated dEs and 
eUCDs suggest that they are a transitional phase between dE,Ns and 
UCDs, with tidal features providing compelling evidence that tidal 
stripping is the linking physical mechanism. In Fig. 4, we show what 
an evolutionary sequence due to tidal stripping would look like in 
observations. The sequence provides a panorama of dE,Ns at different 
stages of such transformation, from anormal dE,N to a UDG, then an 
eUCD (for example, NGVS-UCD509) and, finally, to a ‘naked’ UCD/GC 
(for example, NGVS-UCDS07). 

Ahistory of recent interactions for strongly nucleated dEs and eUCDs 
is also evident from their strikingly different spatial distributions com- 
pared with ‘naked’ UCDs around M87 (Fig. 1b). The distribution of 
strongly nucleated dEs and eUCDsis less concentrated and much more 
inhomogeneous than the bulk of confirmed single-component UCDs 
and GCs”. Strongly nucleated dEs and eUCDs both trace shell-like pat- 
terns and aligned linear substructures around M87 (also see Extended 
Data Fig. 6), mostly in the outer regions of the galaxy, similar to previous 
numerical simulations”. Considering the typical formation timescale 
of eUCDs in the dense cluster environment’, this spatial distribution 
suggests that a group of galaxies has been accreted and fallen into the 
cluster centre in the past 2-3 Gyr (refs. 43,44). 

The nucleation fraction of dwarf galaxies with M, between approxi- 
mately 108-10" M, is higher than 70% (ref. 17), so UCDs represent the 
majority of disrupted normal nucleated dwarf galaxies in this mass 
range, pointing to a steeper galaxy luminosity function in the past®. 
These UCDs may also host a significant number of massive black holes 
in our local Universe*®. Finally, disrupted dE,Ns contribute to the intra- 
cluster light” and stellar halos of massive galaxies. We estimate the 
stellar mass contribution from 16 strongly nucleated dEs, 33 eUCDs and 
89 confirmed UCDs in 1° around M87. Assuming all of these objects are 
stripped nuclei whose progenitors had (fiu) = 1%, they will contribute 
approximately 7 x 10° M,, whichis about 30% of the stellar mass within 
the M87 stellar halo outside of 20 kpc (refs. 45,47). Given that evolved 
tidal streams are extremely faint (usually fainter than 30 mag arcsec *), 
we expect that future deep imaging from Rubin Observatory Legacy 
Survey of Space and Time and the Nancy Grace Roman Space Telescope 
will show more tidal streams around these transitional objects and fully 
illuminate their evolutionary connection. 
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Methods 


NGVS 

The NGVS* is a deep, high-resolution, multiwavelength (u’, 9”, r’, i’ z’) 
imaging survey of the entire Virgo cluster conducted with MegaCam 
on Canada-France-Hawaii Telescope (CFHT/MegaCam”). The survey 
footprint covers 104° (ref. 2) within the virial radii of the Virgo A and 
Virgo B subclusters, centred on M87 and M49, respectively. In the 
g’-band, the point-source detection limit is 25.9 AB mag (100) and 
the extended-source detection limit is 29.0 AB mag arcsec * (2g above 
the sky). Throughout the Article, we assume our objects locate at the 
Virgo distance of 16.5 Mpc, unless the host galaxies can be clearly deter- 
mined and the distance of the hosts are independently known®. 

Taking advantage of the excellent image quality (i’-band seeing less 
than 0.6”, and about 0.8”-1.0” for other bands), Liu et al.* presented 
a sample of 612 confirmed and candidate UCDs in the Virgo Cluster 
using a combination of half-light radius (11 < r, < 100 pc), luminos- 
ity (g’ < 21.5 mag), morphology, colours and available radial velocity 
information. Apart from the UCDs that were already spectroscopically 
confirmed (UCDs with radial velocity v, < 3,500 kms, about one-third 
ofthe sample), an additional 88 UCD candidates in the catalogue were 
selected on the basis of deep u‘g’i’z’K, data to have a high probability 
of being Virgo members. The background contaminant fraction is 
expected to be higher in the outskirts of the Virgo cluster, especially 
in the northernmost region (6 > 16°), where only u’g’z’ colour selec- 
tion is available. Visual inspection by Liu et al.‘ identified 41 objects 
with potential envelopes, including 29 confirmed UCDs and 12 UCD 
candidates, although no structural fitting or photometry of these enve- 
lopes were performed. Strongly nucleated dEs that could potentially 
be UCD progenitors are selected in the NGVS galaxy catalogue“, which 
lists a total of 3,689 galaxies in the Virgo cluster. Note that none of 
the objects in the UCD samples, including eUCDs, are included in the 
dE,N sample derived from the full NGVS galaxy catalogue. VUCD7, 
which is clearly brighter than most UCDs, is included in the strongly 
nucleated dE sample (NGVS 2275) as it was identified as a dwarf gal- 
axy in the NGVS, being quite extended in appearance. Among the 586 
early-type dwarf galaxies that have a total magnitude of 15 < g’ < 21.5 
(-16 <M, < -9.5) and NSC R. < 0.5 arcsec (40 pc, to remove spurious 
nuclei), 51 dE,Ns have f uc > 8%. Galaxies that are identified as compact 
elliptical galaxies (cEs)** were not included in this study. Although 
NSCs are prevalent in nearby galaxies of any morphological type”, 
and stripped nuclear star clusters may come from spiral galaxies, no 
strongly nucleated spiral galaxies were found in Virgo. We expect that 
processes suchas tidal interaction and ram-pressure stripping may have 
transformed those late-type galaxies to early-type galaxies during the 
formation of UCDs, thus such strongly nucleated spiral galaxies may 
be rarein dense environments. We focus on the dwarf elliptical galax- 
ies as they could better show the evolutionary sequence and compare 
their properties with our UCD sample, suchas size, surface brightness 
and colours. 

While the NGVS data has both excellent image quality and deep 
surface brightness sensitivity, we are still limited by the data in both 
respects. The ground-based spatial resolution of CFHT limited the 
parent sample of 612 UCDs to a size cut of r, > 11 pc. This cut removes 
nearly all foreground stars, but also removes UCDs that have smaller 
sizes (as well as many globular clusters). Although more typical 
nuclear star clusters are known to have half-light radius less than 10 pc 
(ref. 17), this falls below the limit for reliable size measurements at 
Virgo distance, which causes the artificial size discontinuity between 
UCDs and normal GCs seen in Fig. 2. Fortunately, UCDs with envelopes 
can have single-component King model r, > 11 pc, even if the decom- 
posed nucleus is smaller. In that way, our data is able to probe to smaller 
nuclear size scales, even if the NSC itselfis unresolved (but the envelope 
is resolved). The NGVS data are thus best suited to finding objects that 
are larger, and have brighter and larger envelopes, that is, those inthe 


earlier phases of tidal stripping. However, it is difficult to decompose 
UCDs that have large (r, > 20-30 pc) and bright NSCs with compact 
envelopes (for example, M60-UCDI1) using only NGVS data. This work 
therefore presents a lower limit on the number of visible stripped nuclei 
inthe Virgo cluster, with fainter and smaller stripped nuclei remaining 
to be found with space telescopes. As a hint towards what remains to be 
found, eight bright globular clusters withg’ < 21 mag also show diffuse 
envelopes (ref. 4, table A1), but are not included in the catalogue of 612 
Virgo UCDs owing to their slightly smaller sizes or mismatched fitting 
results between g’ and i’ band. 


BSDVS 

The BSDVS was conducted with the Case Western Reserve University’s 
Burrell Schmidt telescope”*. The final survey footprint consists of 
15.3 deg” and 16.7 deg’ in the Johnson B- and V-bands, respectively, with 
a pixel scale of 1.45” pixel. The survey covers the main mass concentra- 
tions (for example, the M87, M49, M86 and M84 regions) of Virgo, and is 
optimal for detecting very faint and extended features, including tidal 
features and intracluster lights. The 30 limiting surface brightnesses 
are Up im = 29.5 Mag arcsec “and Hy jim = 28.5 mag arcsec”. 


ACSVCS 

The ACSVCS (GO-9401)* is a program to image 100 early-type galaxies 
of the Virgo Cluster using the Advanced Camera for Surveys onthe HST 
in two widely separated bandpasses (F475W/g’ and F850LP/z’). There 
are 22 UCDs observed by ACSVCS with 0.1” resolution ( full width at 
half maximum), all within approximately 2 arcmin of their host galaxy 
centres™. 


Tidal features 
Various tidal features are commonly found in galaxy cluster environ- 
ments, but here we focus primarily on tidal features that reflect the 
disruption of nucleated dwarf galaxies and the formation of UCDs. 
We performed extensive visual inspections for tidal features around 
all strongly nucleated dEs and all Virgo UCDs on 4 x 4 binned and 
smoothed g’ + r’ + i’ NGVS imaging, and compared NGVS imaging with 
available BSDVS imaging, where there was BSDVS coverage. In addition 
to the previously reported VLSB-A” and VLSB-D”, VCC 1672 and NGVS 
2078 stand out with elongated S-shaped, low-surface-brightness fea- 
tures in their outskirts. Visual inspection of all Virgo UCDs* revealed 
that NGVS-UCD156 is the only UCD associated with a prominent 
20-kpc-long tidal stream that can be seen in smoothed NGVS imag- 
ing. NGVS-UCD156 is in the NGC 4365 group, which is thought to be 
an independent group of galaxies at approximately 6 Mpc behind the 
Virgo Cluster. NGVS-UCD509 near M87 is probably associated with a 
70-kpc-long, but extremely faint, bent tidal stream that is more visible 
onBSDVS imaging. The brightest southern part of this stream, located 
48” (approximately 38 kpc) southwest of NGVS-UCD509, is catalogued 
as NGVSUDG-17 (ref. 30) and may be stripped material from the pro- 
genitor dE,N, which can be seen in Fig. 4. Both NGVS and BSDVS reveal 
the faint but noticeable over 20-kpc-long tidal stream that originates 
froma bright (M, = -11.34 mag) globular cluster NGVSJ1229 + 1343 with 
a diffuse envelope, suggesting NGVSJ1229 + 1343 is a remnant nuclear 
star cluster of a disrupting dwarf galaxy. Located near Markarian’s 
Chain, NGVSJ 1229 + 1343 is at similar projected distances to NGC 4459 
(16 arcmin = 75 kpc), NGC 4477 (19 arcmin = 90 kpc) and NGC 4473 
(23 arcmin = 110 kpc), so it is difficult to unambiguously identify an 
obvious host that may be causing the disruption of NGVSJ1229 + 1343. 
Another Virgo nucleated UDG, NGVSUDG-O1 (ref. 30), is embedded in 
avery long, narrow filamentary tidal stream” that provides additional 
evidence for the UDG-UCD formation connection. The tidal features 
around one bright strongly nucleated dE, VCC 1148 (NGVS 1781), and 
its companion dwarf galaxy VCC 1153, have also been reported”. 
Most of the remaining strongly nucleated dEs (47 of 51) donot show 
signs of tidal features down to an imaging depth of uy ~ 29 mag arcsec’, 
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but this does not preclude there having been a history of interactions. 
Anorbital study” of one nucleated UDG (VCC 615) without tidal features 
suggests that it likely passed near the cluster centre and thus experi- 
enced some degree of tidal stripping, but the resulting tidal features 
may be too faint to be detected with current NGVS and BSDVS imaging. 
The elongated GC distribution and dynamical mass measurement for 
some UDGs are also independent evidence for the tidal stripping”, 
although we cannot rule out the possibility that some strongly nucle- 
ated dEs directly formed an over-massive NSC, especially at the faint 
end, where f,,,. is generally higher than for dE,Ns at higher masses. 
However, the continuum seen in fuc and effective radius, and the very 
different spatial distribution for strongly nucleated dEs compared to 
more typical dE,Ns (Extended Data Fig. 6), suggests that these galaxies 
are part of an evolutionary sequence. Meanwhile, we noticed that many 
strongly nucleated dEs in the Virgo core regions lack GCs, whichis also 
consistent with the tidal stripping scenario. If they were not formed 
via stripping then an explanation would be needed for the formation 
of such a massive central NSC, one containing more than 10-40% ofa 
galaxy’s stellar mass. 


Two-dimensional image decomposition and identification of 
eUCDs 

We used GALFIT®” to fit multiple models to each UCD in 200 x 200 
pixel (0.62 x 0.62 arcmin) subimages: a single King® (K), a single Sérsic™* 
(S), and two-component King + Sérsic (K + S), double Sérsic (S + S) or 
point spread function (PSF) + Sérsic (P + S) models. The effective PSF 
was generated by the Astropy package Photutils® using several nearby 
unsaturated and isolated stars. For objects that had both NGVS and HST 
imaging, we found that for the two-component models, the K + S were 
more consistent between the two datasets than S + S. The P + S models 
often had worse results, even in NGVS imaging, indicating that many 
NSCs are slightly resolved. For these reasons, we subsequently only 
discuss the K + S two-component fits. 

The initial fitting parameters are based on the magnitude and 
half-light radius from previous measurements‘. For the initial Sérsic 
index (n) and King concentration (c = log,o(r,/r,), where r, is the tidal 
radius and r, is the core radius), we use typical values of n = 2 and 
c=1.5, which are broadly representative of early-type galaxies and 
globular clusters, respectively. For the two-component fit, we adopt an 
initial guess of f,,,.=30%, and outer Sérsic R, = 100 pc and n =1.5. The 
input shape of all components is set to an axis ratio q = 1, given the fact 
that most UCDs are fairly round. All model parameters are allowed to 
vary in the fitting, except a fixed shape parameter a = 2 for the King 
model, plus a constraint of Sérsic index between n = 0.5 and n = 8.0. 
The x and y position are fixed to be the same for two-component fit- 
ting. The NGVS images are already globally sky-subtracted, and any 
residual local sky background and gradient was fitted and subtracted 
with GALFIT before fitting the UCD, whichis performed after masking 
all objects in the image. When the UCDs are severely influenced by 
the non-uniform stellar light from nearby galaxies (about 20 UCDs 
in our 612 UCDs have this problem, but most of them are covered in 
ACSVCS), the neighbouring galaxies were first modelled and sub- 
tracted, although the higher noise related to the bright sky back- 
ground makes it more difficult to identify stellar envelopes in these 
objects. Generally, the single-component fit converged quickly, and 
residual images show that UCDs are well modelled and subtracted. 
The typical Sérsic index (n) for a single-component fit is n = 1-4, 
with a median of n = 2.58. The typical King concentration index fora 
single-component fit is c = 1-2, with a median of c = 1.54. A small but 
important fraction (approximately 15%) of objects, however, cannot 
be reproduced well with a single profile. Under these circumstances, 
the Sérsic index usually reaches our upper limit of n = 8, which implies 
that these objects may have extended envelopes. Objects that have 
a large Sérsic index in the single Sérsic model also show unphysical 
parameters in their single King model, including a runaway large tidal 


radius (r,), large concentration index c (for example, c = 3), and very 
high central surface brightness. In these cases, two-component models 
provide a more physically meaningful description of their surface 
brightness profiles. 

In this work, we address a potential structural continuum from 
objects with two clearly distinct components to those with a single 
component, and are always limited by the resolution of ground-based 
imaging. Although there is no fully objective or completely automated 
procedure to determine the existence of envelopes and select the 
best-fit model, we use a uniform methodology that is based on exam- 
ining all the information from the data and models, plus our reasonable 
expectations for the structures of UCDs, case-by-case. Here we provide 
a summary of our classification of eUCDs. Our initial judgement is 
based on the single-component King and Sérsic model fits. We choose 
these models because, unlike, for example, power-law profiles, they 
are physically motivated, always integrate to a finite luminosity and 
are expected to fit a wide range of stellar dynamical systems. The sin- 
gle Sérsic index and single King profile concentration index served as 
first selection criteria for eUCDs. In practice, most typical eUCDs have 
Sérsic index n > 6 and c > 2 through our fitting. For these objects, we 
performed the two-component fit by randomly generating reason- 
able initial guess parameters many times for each model choice and 
fine tuning the parameters based on its one-dimensional (1D) surface 
brightness profile until we got the best-fit result. The goodnesses of 
thetwo-componentfits are examined by comparing the reduced x’ and 
residual maps with the single-component fits. A good two-component 
fit is expected to subtract all the extended light excess, whereas a 
single-component fit will still leave systematic residual features (for 
example, outer rings). We further check the structural parameters of the 
best two-component fit result. These parameters should be reasonable 
and consistent with the general properties of dE,Ns and UCD envelopes. 
We exclude objects if the component axis ratio q = b/a < 0.3, or if the 
outer component has a higher central surface brightness than the inner 
component. For the rest of the UCDs, for those that have a more realis- 
tic single Sérsic index and single King profile concentration index, we 
visually inspect their 1D surface brightness profile and two-dimensional 
(2D) images, and search for any noticeable inflections in the profile 
that are not caused by nearby sources. We then repeated the selection 
process above. The relative goodness between the one-component 
and two-component fit, and the choice of the best-fit model, was also 
assessed using the Bayesian Information Criterion (BIC®). The inte- 
grated colour information of the inner and outer components (shownin 
Fig. 3 and Extended Data Fig. 7) was derived by fixing the best-fit g’-band 
structural parameters and fitting only the magnitude in other bands. 
Examples of fitting results for eUCDs in Fig. 1 are shown in Extended 
Data Fig. 2. 

All UCDs and UCD candidates are classified into four categories 
(0, no envelope; 1, possible envelope; 2, probable envelope; 3, certain 
envelope) according to the prominence of stellar envelopes and fit- 
ting result, which are cross-checked by experienced team members. 
We found the trend in median difference of BIC scores between the 
one-component and two-component fits (ABIC = BIC,,,. ~ BIC,,,,) is 
consistent with our classification, for example, ‘certain’ eUCDs tend 
to have the highest positive ABIC, while ABIC of ‘possible’ eUCDs are 
smaller and are more consistent with UCDs with no envelopes. We 
selected 84 objects to have at least ‘possible’ envelopes (categories 1,2 
and 3), including 45 that were already confirmed to be Virgo members 
with radial velocity measurements from previous studies”*”’ and 39 
UCD candidates that needed spectroscopic follow-up. We observed 
38 of these 39 with Gemini North/GMOS in 2021 and 2022, and we 
confirmed 17 bona fide Virgo UCDs (see ‘Spectroscopy’ section). At 
this point, we have 50 ‘certain’ and ‘probable’ eUCDs, and 12 out of 
13 ‘possible’ eUCDs, that are confirmed Virgo members, eliminat- 
ing nearly all background galaxy contaminants. ‘Possible’ eUCDs are 
typically UCDs that have a high Sérsic index or King concentration 


but for which a stable two-component decomposition is difficult to 
obtain. In these cases, both the structural parameters and colours 
are expected to have significant uncertainties due to model degen- 
eracy. For simplicity and consistency, we only plot the results of 
‘certain’ and ‘probable’ eUCDs as eUCDs, and treat ‘possible’ eUCDs 
as one-component UCDs throughout the Article. Not including the 
13 ‘possible’ eUCDs would not change any conclusions of our Article 
(spatial distribution, evolutionary sequence and so on). ‘Possible’ 
eUCDs may represent more highly stripped nuclei that are evolving 
to ‘naked’ UCDs. Finally, we analysed 22 UCDs with HST imaging in 
the ACSVCS (see ‘Consistency Checks’ below). Of these, five objects 
that were classified as ‘No envelope’ using NGVS imaging were subse- 
quently found to have envelopes when analysed with the benefit of HST 
resolution. These five additional eUCDs already had measured veloci- 
ties and were confirmed Virgo members, bringing our total sample 
of eUCDs to 55. 


Uncertainties on derived photometric parameters 

The main uncertainties in our eUCD decomposition include model 
mismatch, sky background subtraction and PSF variation. Previous 
work has shown that GALFIT usually underestimates errors”. To get 
more realistic measurement errors, for each eUCD we generate 100 
realizations of mock eUCD images based on derived best-fit param- 
eters accounting for the Poisson noise associated with the source, 
background noise and read noise, which are calculated based on the 
properties in specific images. We also take into account the possibil- 
ity of asystematic sky subtraction error by introducing a sky offset 
sampled from a Gaussian with a sigma of 3%, a value that probably 
overestimates the sky subtraction bias under our detailed analysis 
but is used to be conservative. We then use GALFIT to repeat the 
decomposition for these mock UCDs. The 1o uncertainties for size, 
magnitude and colour are calculated using the distribution of the 
fitted parameters for these mock observations. For the magnitude 
and colour uncertainties we add in quadrature the derived systematic 
uncertainties for NGVS photometry“. A typical error bar for size and 
luminosity is shown in the Fig. 2 legend. The typical magnitude error is 
around 0.15 mag, and the envelope size measurement error is around 
15%. Colour and magnitude uncertainties are shown for all strongly 
nucleated dwarfs and eUCDs in Fig. 3. Although the uncertainties are 
unique to each object and may be influenced by nearby neighbours or 
varying backgrounds, the typical uncertainties in both size and colour 
are smaller than the scatter seen in the sample. We notice the colour 
uncertainties of bright and compact UCD envelopes are less affected 
by sky subtraction error than those much more diffuse strongly nucle- 
ated dEs, for which the sky subtraction error is considered to be the 
dominant uncertainty. 


Structural decomposition of tidally disrupting dE,Ns 

We ran GALFIT on these galaxies to derive the structural parameters 
of the objects ing’-band. We performed accurate local sky subtraction 
using the GNU astronomy utility NoiseChisel®, which detects extremely 
faint structures, allowing us to create masks of low-surface brightness 
features in the image, including nearby galaxies and optical ghosts 
around bright stars, as well as artefacts like diffraction spikes. We then 
estimated the 2D local background map from unmasked regions using 
Photutils Background2D. 

Allnearby objects, including background sources overlapping with 
the target galaxies, are detected and masked using segmentation maps 
generated by Source Extractor”. Masks were created in two stages. 
The primary masks were created based on the original images. We 
fed the primary images into GALFIT to model and subtract the target 
galaxy. In the second stage, the residual map of the initial fit was used 
to create a second mask. The final mask for a given image is the union 
of these two masks. This process is done independently in u’-, g’- and 
i’-band images, and the final image mask is the union of the masks in 


all three filters. Any bright objects other than the smooth component 
and central nucleus of the galaxy are masked in the final image. 

We took an iterative approach to our structural analysis. Instead 
of fitting several components simultaneously, we first fit the outer, 
extended component on co-added images. This component is most 
sensitive to the accuracy of the background subtraction and the Sérsic 
index, so we initially fix the Sérsic index to n = 1 and fit for R, and the 
total magnitude of this component. All foreground and background 
objects, including nuclear components, were masked in this fit. We then 
performed an unconstrained second fit on the residual images, focus- 
ing onthe profile of inner diffuse envelopes and nuclear components. 
Lastly, we put the parameters from extended components and nuclear 
components together and re-ran the fit simultaneously on images in 
multiple bands (allowing all parameters to be unconstrained) until 
we arrived at the best description of their overall profiles. The best-fit 
result for disrupting dE VCC 1672 can be found in Extended Data Fig. 2. 

Using three components to decompose the structure of tidally 
disrupting dEs rather than two components is motivated not only by 
the data but also by theoretical and numerical studies that have sug- 
gested a more complex structure during tidal stripping”*””. Part of 
the stripped, unbounded stars will relax to virial equilibrium after the 
most recent pericentric passage, which causes a ‘tidal break’® that cor- 
responds to the time elapsed since pericentre. Actually, simulations of 
UCD formation via tidal stripping also present a similar ‘tidal break’ in 
the surface brightness profile during the early stages”. 

Moreover, previous observations” of disrupting dE,Ns have already 
shown similar NSC + envelope + tidal stream structures. More com- 
plex structures are expected to be generated by multiple pericentric 
passages, which may help us to constrain the orbital history of these 
disrupting systems, for example, VLSB-A shows a potential NSC + 
envelope + elongated bar-like structure + ultra-diffuse outer compo- 
nent. The stripped stars would eventually disperse with time, and the 
outermost component gradually becomes invisible. The net effect is 
consistent with the evolutionary sequence we observed in Fig. 4. 


Spectroscopy 

We obtained radial velocities for 9 dE,N candidates and 38 eUCD candi- 
dates using the Gemini Multi-Object Spectrograph” on the 8 m Gemini 
North telescope during the 2021 A and 2022 A semesters. We used the 
0.75 arcsec long-slit and B600 grating for 2021 A, with a central wave- 
length at 5,250 A and a spectral resolution of R = 1,125. For 2022 A, we 
used the 1.00 arcsec long-slit, witha central wavelength at 6,250 A. For 
g =21.5 mag objects, the signal-to-noise (S/N) per resolution element 
is approximately 5 with 45 min of exposure time. 

We reduced the Gemini Multi-Object Spectrograph spectra using 
the Pypelt® package. After extracting the 1D spectra, we derived radial 
velocities with penalized pixel fitting software pPXF”””. We used high 
signal-to-noise (S/N = 80) MMT/Hectospec stacked spectra of blue 
and red Virgo GCs as velocity templates”. We used the wavelength 
ranged = 5,500-6,800 A, containing the strong H, (6,563 A) absorption 
feature, whichis among the most prominent spectral features in the 
relatively low S/N spectra we obtained. We identified 17 eUCDs as Virgo 
members with radial velocities V, < 3,500 kms ‘and velocity uncer- 
tainties Av = 20-90 km s”, anda median Av = 35 km s”. The remaining 
21 eUCD candidates are background galaxies, with redshifts ranging 
fromz= 0.028 to 0.23. The high contamination fraction shows the key 
role of the spectroscopic survey in making an unbiased interpreta- 
tion of properties of eUCDs. The bona fide eUCDs and background 
galaxies confirmed in the spectroscopic survey show obvious dif- 
ferences in spatial distribution and shape. Background galaxies are 
mostly found in the low-density outskirt regions of Virgo, and alleUCD 
candidates with elongated (b/a < 0.6) shapes were found to be back- 
ground galaxies. All confirmed eUCDs are round with (b/a 2 0.7). All 
nine observed dE,Ns are Virgo members, consistent with the fact that 
their surface brightnesses and sizes are similar to other Virgo galaxies. 
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We expect that the contamination fraction for the rest of the strongly 
nucleated dEs without spectroscopy is also low. Example spectra for 
a bona fide eUCD and a background galaxy are shown in Extended 
Data Fig. 3. 


Consistency checks 

We fitted 22 UCDs that are also covered in the ACSVCS® and com- 
pared the results with NGVS fitting. Empirical PSFs in the F475W and 
F850LP filters were derived using DAOPHOT II (ref. 78) using archival 
observations of fields in the outskirts of the Galactic globular cluster 
NGC 104 (47 Tucanae)”. For galaxies that have good sky subtraction in 
both NGVS and ACSVCS, we find their best-fit model and parameters 
are generally consistent, including effective radius, Sérsic index and 
King concentration. For example, NGVS-UCD190, NGVS-UCD330 and 
NGVS-UCD395 are three eUCDs identified in NGVS imaging, which can 
be more easily decomposed into two components in HST imaging. The 
fauc Values derived from NGVS imaging are generally about 5% lower than 
those from HST imaging. The envelope component when fitting NGVS 
imaging tends to havea higher Sérsic index, and the effective radius is 
10-20% larger than when fitting HST images. These behaviours can be 
explained by the shallower depth of HST imaging, which misses a small 
fraction of light in the outer envelope with surface brightness between 
26 and 28 mag arcsec”. This results in a smaller effective radius, lower 
Sérsic index and a slightly higher f u estimation when using HST imag- 
ing. A comparison between HST imaging and NGVS image fitting for 
an eUCD is shown in Extended Data Fig. 4. 

In addition to the eUCDs mentioned above that were in the 
ACSVCS, but could be identified in the NGVS as eUCDs, we were also 
able to use HST resolution to identify additional eUCDs with com- 
pact envelopes that were unresolved in the NGVS. Five UCDs in the 
ACSVCS have a single-component Sérsic index n > 4, in which three 
have n>7 (NGVS-UCD167, NGVS-UCD386 and NGVS-UCD822). 
The two-component fits to these objects suggest that NGVS-UCD167 
has aR, = 30 pc envelope and NGVS-UCD822 has aR, = 24 pc envelope, 
whereas NGVS-UCD386 shows a power-law profile and we are unable to 
decompose this object. NGVS-UCD192 has a prominent R, ~ 37 pc enve- 
lope, and NGVS-UCD548 also show a two-component profile, with enve- 
lope size R, = 35 pc. Another eUCD from HST imaging only is M60-UCD1 
(NGVS-UCD753), which has a compact envelope (R, = 50 pc), consist- 
ent with previous measurements”. Our NGVS work finds eUCDs with 
envelopes having R. > 40 pc, whichis probably a selection limit based 
on spatial resolution. These eUCDs with envelopes having R, < 40 pc 
can only be resolved with HST imaging. A comparison between HST 
imaging and NGVS image fitting for an eUCD discovered only in HST 
imaging is shown in Extended Data Fig. 5. The eUCD-to-UCD fraction 
in NGVS imaging is over 12%. We suggest that this is a lower limit, 
and that the actual fraction should be higher. However, the current 
ACSVCS sample may be biased in that the ACSVCS was a targeted gal- 
axy survey, and so preferentially imaged areas where UCDs are more 
likely to form through interaction with a larger host. Future wide-area 
surveys of the Virgo cluster with HST, Euclid or the Roman Space 
Telescope will be needed to build an unbiased sample of UCDs/GCs 
with smaller envelopes. 


Data availability 


NGVS data can be accessed via the Canadian Astronomy Data Centre 
(CADC). ACSVCS data is available via the Mikulski Archive for Space 
Telescopes (MAST) with the programme ID HST GO-9401. Data products 
of Burrell Schmidt Deep Virgo Survey can be downloaded from http:// 
astroweb.case.edu/VirgoSurvey/. Gemini/GMOS spectroscopic data of 
eUCD candidates taken during 2021-2022 can be downloaded through 
the Gemini Observatory Archive with programme IDs GN-2021A-Q-208, 
GN-2022A-Q-206 and GN-2022A-Q-307. Source data are provided with 
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Extended Data Fig. 1|2D decomposition of eUCDs. First column: Radial 
surface brightness profile of the UCD and best-fit result of the two-component 
fitting. Alle UCDs can be well-fit with an inner King model and an outer Sérsic 
model. The total, King, and Sérsic model profiles are shown with red, magenta 
and green dashed lines, respectively. The grey dashed line is the PSF profile. 


Error bars represent louncertainties. The fitting residuals are shown below 

the source profile. Second column: NGVS u‘g’i’ colour composite image. 

Third column: Theg’-band image of UCD, with contours showing the surface 
brightness at 4, = 24 to 27 mag arcsec”. Last column: (data-~model) residual map. 
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Extended Data Fig. 2| The best-fit three-component model of VCC 1672. 
The right panels display, from top to bottom, images of the data overlapped with 
isophotes, the best-fit model, and the residuals. North is up and east to the left. 
The left panels display the isophotal analysis of the 2D images and model fitting. 
From top to bottom, the panels show radial profiles of the axis ratio (q), the 
position angle (PA), the g’-band surface brightness, and the fitting residuals. 
For an illustration of the breaks in position angle and axis ratio, the effective 
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radius of envelope components and outer components are plotted as vertical 
dot-dashed lines. Ellipses represent the effective radius of each individual 
component shown in the model image, which uses the same colour legend asin 
the surface brightness plot. We note that the shifts in both q and PA roughly 
correspond to the radius at which the outer component begins to dominate 
over the central envelope, at Re ~ 500 pc. Error bars represent louncertainties. 
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Extended Data Fig. 3 | Examples of Gemini/GMOS spectra of a bona fide objects, both of which show extended envelopes. The strong H, absorption is 
eUCD anda background galaxy. The Gemini/GMOS optical spectrum (black), the main diagnostic for bona fide eUCDs. In contrast, the measured redshift of 
the best-fitting template (red) and fitting residuals (gray) are plotted in the left NGVS-UCD147 is z ~ 0.133, a background galaxy composed of bulge+disk. 
panels. Right panels show the g’-band NGVS imaging and contours for each 


NGVS-UCD3g0. 7 $ Data 


9 


-- King+Sérsic 


100 
Radius (pc) 

= > Data 
N i 2. 
6) -- Kingt+Sersic 
a PSF 

1S) 

= 

© 
~ 

od) 

© 

= 

= 

in 

~ 

+ 

LL 

3 
= 
in 
a oo 0290 
ire ©0000000000009°05,0009 © 

a. Sočo 
10 100 
Radius (pc) 
Extended Data Fig. 4 | A comparison between NGVS and HST imaging for shows marginal two components in NGVS imaging, and can be much clearly 


NGVS-UCD330. The upper and lower row are NGVS and HST data, respectively. decomposed into two components in HST imaging. The total, King and Sérsic 
Contours in the left column represent the surface brightness at u,,.=24-27mag model profile are shown with red, magenta and green dashed line. The grey 
arcsec? in NGVS and H azsw = 23 - 25 mag arcsec” for HST. In the right column dashed line is the PSF profile. Error bars represent louncertainties. 

is the one-dimensional surface brightness profile of the UCD. NGVS-UCD330 
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Extended Data Fig. 5| A comparison between NGVS and HST imaging for NGVS-UCD167. NGVS-UCD167 shows only one component in NGVS imaging, but can 
be resolved to two components in HST imaging. 
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Extended Data Fig. 6 | Spatial distribution of compact stellar systems 

in the Virgo cluster. Yellow and dark red circles represent UCDs and eUCDs, 
respectively. Blue circles and squares are strongly-nucleated dEs with and 
without velocity measurements. Gray circles represent all UCD candidates 
without velocity measurements‘. Disrupting dEs”"***° and eUCDs that show 
tidal features are labelled with green stars. The Gaussian-smoothed distribution 
of Virgo dE,Ns identified in NGVS is shown in the background. While the density 


peak ofall dE,Ns is located around M84 & M86 (northwest of M87), most eUCDs 
and strongly-nucleated dEs are concentrated around M87. Massive Virgo 
galaxies (M,.<-20) that could be primary contributors of tidal disruption are 
shown in white dashed circles, with circle sizes corresponding to their total 
luminosity. Most strongly-nucleated dEs and eUCDs have clear associations 
with larger galaxies, or form coherent substructures. 
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Extended Data Fig. 7 | Nuclear star cluster colour-magnitude diagram for colour-magnitude relation, the nuclear star clusters in UCDs have colours that 


are very similar to UCDs without envelopes, and the NSCs of normal dE,Ns. 10 


eUCDs, strongly-nucleated dEs, and Virgo early-type galaxies. Unlike the 
error bar for the colour and magnitude is shown for each object. 


colours of the stellar envelopes, which are offset from the main galaxy 


